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Abstract. The observed correlations, between the characteristic ages and dipole sur- 
face magnetic field strengths of all pulsars, can be well explained by magnetic field 
decay with core temperatures of 2 x 10* K, ~ 2 x 10^ K, and ~ 10^ K, for magne- 
tars, normal radio pulsars, and millisecond pulsars, respectively; assuming that their 
characteristic ages are about two orders of magnitude larger than their true ages, the re- 
quired core temperatures may be reduced by about a factor of 10. The magnetic decay 
follows a power-law and is dominated by the solenoidal component of the ambipolar 
diffusion mode. In this model, all NSs are assumed to have the same initial magnetic 
field strength, but different core temperature which do not change as the magnetic field 
decays. This suggests that the key distinguishing property between magnetars and nor- 
mal pulsars is that magnetars were born much hotter than normal pulsars, and thus have 
much longer magnetic field decay time scales, resulting in higher surface magnetic field 
strength even with the same ages of normal pulsars. The above conclusion agrees well 
with the observed correlations between the surface temperatures of magnetars and other 
young NSs, which do not agree with the cooling dominated evolution of neutron stars. 
This suggests a possible scenario that heating, perhaps due to magnetic field decay, 
balances neutron star cooling for observed pulsars. 



1. Introduction 

The strength of the surface dipole magnetic field of a neutron star may be estimated 
directly from the period and period derivative, i.e. B^ w 3.3 x \Qi^^{PP)^l^ G, which 
assumes that its rotational energy is lost entirely via magnetic dipole radiation. Then the 
spin-down age (characteristic age) can also be estimated as Tc = PI2P. Then all pulsars 
can also be displayed in the - B diagram (Figure 1), instead of the P - P diagram. 
Clearly all pulsars can be divided into three clusters in Figure 1 ; from top to bottom, 
the three cluster are magnetars, normal pulsars, and millisecond pulsars, respectively. 

Magnetars, which are normally referred to as neutron stars with surface magnetic 
fields larger than the quantum critical value, Bcr - m^c^ /eh « 4.4 x 10^^ G, have been 
proposed as the most promising energy source for Anomalous X-ray Pulsars (AXPs) 
and the Soft Gamma-ray Repeaters (SGRs). In the - B diagram, magnetars are 
composed of a distinctive group with higher magnetic field (~ 5 x 10'^ - 2 x 10^^ G) 
and younger characteristic ages (~ 10^ - 10^ yrs). The two classes of X-ray Pulsars, 
AXPs and SGRs, differ from the common accretion-powered pulsars in massive X- 
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Figure 1. The - B diagram. Blue triangles: SGRs and AXPs. Magenta 
circles: other pulsars, including normal radio pulsars, radio-quiet pulsars, high 
energy pulsars, millisecond pulsars, binary pulsars, etc. Thin lines: the best 
fitted line for magnetars (dotted line), ordinary pulsars with dipole field strength 
4 X 10'° - 2 X 10'^ G (dashed line), and millisecond pulsars with dipole field strength 
< 10^ G (dot-dashed line). Thick solid lines: 'death line' (black), magnetic field de- 
cay for = 2 (blue), - 0.2 (red), T% - 0.001 (purple). For magnetars, their char- 
acteristic ages (Tc) and surface magnetic fields Bjurf (^d) are from McGill SGR/AXP 
Online Catalog (http://www.physics.mcgill.ca/~pulsar/magnetar/main.html); 
for other pulsars, and Bjurf are from the ATNF I'ulsar Catalogue 
(http://www.atnf.csiro.aU/people/pulsar/psrcat/j . 



ray binaries mainly in the following observational properties (see iMereghettil 12008 , 
for a review): they probably have no binary companions, but their persistent X-ray 
luminosity can be larger than their spin-down power; they have periods of activities 
during which they emit numerous short bursts in hard X-ray/soft gamma-ray band, and 
the bursts have peak luminosity up to ~ 10^^ erg s~' and last for ~ 0.01 - Is; they 
have a short hard spike, with luminosity larger than ~ 5 x 10'*^ erg s~\ followed by a 
long pulsating tail, the so-called giant flares, a very striking phenomenon observed only 
from SGRs. 

It is believed that the energy reservoir fueling the SGR/AXP activity is their ex- 
treme magnetic fields ([Thompson & DuncanI 1 19951. IT996. hereafter TD96). For in- 
stance, their magnetic fields have enough free energy to power the giant flares; mean- 
while the short duration of their initial spikes is consistent with the pro pagation with 



Alfve n speed of the magnetic instability over a whole neutron star surface ([Thompson & Duncan 



19951). and the niagne tic confinement of the hot plasma is responsible for the pulsating 



tails (IMereghettil l2008b . However, the recent discovery of a SGR with a very low sur- 
face dipolar magnetic field, SGR 0418-1-5729, whose dipolar magnetic field cannot be 
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greater than 7.5 x 10'^ G, implying that a hig h surface dipola r magnetic field is not nec- 
essarily required for magnetar-like activities ( Rea et all2010f) . And the outer gap model 
of magnetars may also be challenged by the null detection in a Fermi/LAT observation 
of AXP 4U 0142+61 (IXong et al.(l2O10l) . Thus, it is likely that t he magnetar activity is 
driven by the magnetic energy stored in the internal toroidal field ('Thompson & Duncan] 
\V995, 2001) or the surface multipolar field. After the discovery of SGR 0418+5729, 
the magnetar population includes objects with a wider range of B-field strengths and 
ages; however their evolutionary stages still remain far from clear. 

In this report, by adopting the magnetic field decay model to the magnetars, normal 
pulsars, and millisecond pulsars, we show that the inner core temperatures of the three 
types of pulsars are systematically different, consistent with the observed three clusters 
shown in Figure 1 and the observed surface temperature of magnetars and other young 
normal pulsars. It is found that the surface and core temperatures of magnetars are the 
highest (include the low magnetic field magnetar SGR 0418+5729) and remain constant 
for at least 24 Myrs. 



2. Magnetic Field Decay of All Pulsars 



Goldreich & Reiseneggej d 19921. hereafter GR92) studied several avenues for magnetic 



field decay in isolated neutron stars: ohmic decay, ambipolar diffusion, and Hall drift. 
Depending on the str ength of the magn etic field, each of these processes may dominate 
the evolution. Hevl & Kulkamil (I1998'. hereafter HK98) applied the model to magne- 
tars, and found that for sufficiently strong nascent fields, field decay alters significantly 
the cooling evolution relative to similarly magnetized neutron stars with constant fields. 
Following GR92 and HK98, a simple differential equation can be used to describe the 
dipole magnetic field decay (HK98) 



dBp 111 
— ^ = -Bp( + + - 

dt ^ohmic ^ambip ^hall 



-Bp(- + + — ), (1) 



where Bp is the strength of the magnetic field at the pole at the surface of a neutron 
star, and fohmio ^ambip> and fhaii are the decay time scales due to ohmic decay, ambipolai" 
diffusion, and Hall drift, respectively. Taking ctq = 4.2 x 10^^rg"^(p/pnuc)^s~', where 
denotes the temperature in units of 10^ K, pnuc = 2.8 x lO''* g cm~^, we have fohmic ~ 
2 X 10^1 yr, where L5 is a characteristic length scale of the flux loops through 

ig Pnuc 

the outer core in units of 10^ cm. The time scales for Ambipolar diffusion are given by, 
Cbip ~ 3 X 10^:^ yr, and t^^^,^ ~ ^(1 + 5 x lO'^L^r^) yr, where Bn is the field 
strength in units 10'^ G, the f^[,jp represent for the time scale of solenoidal component 
and f!|'i^(,;p for irrotational component; note that the latter is only for the case of modified 

URCA reactions. Finally the time scale of the Hall cascade is fhaii ~ ^^^^^^^(]f~^ 

To compare eq. (1) with data shown in Figure 1, we take L5 = 1, p/pnuc = 1> 
initial field strengths Bp - 10'^ G, and begin the integrations at f = 0.1 yr, and assume 
constant core temperatures for the three classes of pulsars. In Figure 1, we plot the cal- 
culated magnetic field decay histories with eq. (1) f or three typical core tem peratures; 
for comparison the pulsar 'death line' is also plotted (IChen & Rudermanll 19931) . Several 
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Figure 2. Observations of the surface temperatures of magnetars and young NSs. 
Left panel: surface temperature vs. characteristic age. The red soHd line is a ba- 
sic theoretical coohng cur ve of a nonsuperfluid neutron star with M = 1.3 Mq 
(lYakovlev & PethicS 12004 *). Right panel: surface temperature vs. magnetic field. 
The green dotted line and dashed line are heating balance lines (HBL) for different 
parameters, and the band between them is heating balance area. The red solid line 
is a theoretical curve for magnetothermal evolution of magnetized neutron stars, in- 
cluding crustal heating by magnetic field decay (APM08). Each number marked in 
the plot is the characteristic age of the pulsar in unit of 10^ yrs. 

observations can be obtained from Figure 1: (i) Magnetars also follow the 'death line' 
(albeit with low statistics), calculated assuming the radio pulsar mechanism, (ii) The in- 
ner core temperature of magnetars is systematically higher than other types of pulsars. 
The core temperature of magnetars has a very narrow range and is about 2 x 10^ K. 
The core temperatures of ordinary pulsars distribute in a relatively wide range, from 
2 X 10^ K to 2 X 10^ K; the center is about 2 x lO' K. The core temperature of miUisec- 
ond pulsars is around 10^ K. (iii) Several radio-quiet pulsars and normal radio pulsars 
fall right on the curve for magnetars, and thus may have the same core temperature. 

All observed NSs appear in the power-law parts of the model curve s, i.e., fisurf 
T~°'^. As we have shown in a companion paper in the same proceedings ("Zhang & Xie| 
201 lb . this corresponds to the solenoidal component of the ambipolar diffusion dom- 
inated decay mode with a constant core temperature. Interestingly, exactly the same 
relation is obtained in the magnetic dipole spin-down model with a constant P, i.e., 
P <^ \. This immediately suggests that the observed extremely small P of NSs is due 
to magnetic field decay with constant core temperature. 



3. Core Temperatures of Magnetars and Young Pulsars 

Currently the core temperature of a NS can only be inferred from its surface temperature 
Tj. The spectra of some NSs are well-described by a simple blackbody (BB) radiation 
which is believed to be its thermal surface emission. For NSs with spectra that include 
two BB components, we adopt the temperature of the component which dominates the 
spectrum. In Figure 2 we show the observed - and T^- B relations for magnetars 
and young NSs. In the left panel of Figure 2, it is shown that the surface temperature 
distribution can be divided into two separate groups. The normal radio pulsars and 
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several radio-quiet pulsars may follow the cooling curve with large scatters, though a 
constant temperature may also describe the data almost equally well. However, the 
surface temperatures of magnetars (> 4 x 10^ K) are much higher than that of normal 
young NSs and also almost constant for at least twenty million years, consistent with 
the model parameters used in Figure 1. 

IPons & Geppert (2007) reported a correlation between the surface temperature and 
the magnetic field and investigated similar diagram as the right panel. They approxi- 
mated the correlation as J^g - CB^ j^, the so-called heating balance line (HBL), (see 
also APM08), where C is a constant that depends on the thickness of the crust with a 
typical value of C ^ 10. They argued that the thermal evolution of NSs with B > 10^^ G 
is predominantly determined by the strength of the magnetic field, and NSs with field of 
~ 10^^ G should cool much more rapidly than those with field of ~ 10^^ G and higher. 
They further suggested that "A NS will begin its life high on the figure with some B^. 
As it cools it moves vertically downward, until decay of its field causes the trajectory to 
bend to the left. The star eventually reaches its HBL, and then continues moving down 
it." 

However, the above scenario is not supported by data; the data are more con- 
sistent with a scenario that both magnetars and normal pulsars move horizontally to 
the left, but with initially very different temperatures, i.e., their magnetic fields de- 
ca y but core t emper atures remain the same, again consistent with the results in Figure 
1. Rea et al. (I2010h suggested that the magnetic energy stored in the internal toroidal 



field that power the violent activities for magnetars, rather than surface dipole field. 
Thompson & DuncanI ( 199dh investigated the process that magnetic field decay can pro- 



vide a significant source of internal heating, and the possibility of obtaining an equilib- 
rium between neutrino cooling and heating through magnetic field decay for B ~ 10'^ G 
and r ~ 10*^ K. The internal toroidal field probably plays this role and therefore keeps 
the core temperature as well as surface temperature of magnetars nearly constant. 



4. Summary and Discussion 

In this report, we found that the observed - B correlations of all pulsars can be well 
explained by magnetic field decay with core temperatures of 2 x 10^ K, ~ 2 x 10' K, 
and ~ 10^ K, for magnetars, normal radio pulsars, and millisecond pulsars, respec- 
tively; the decay is dominated by the solenoidal component of the ambipolar diffusion 

mode with a time scale of f^^^^^ ~ 3 x 10^ -^y^ yr. With just one free parameter (the 

core temperature), all data of pulsars can be reasonably well described. We emphasize 
that in Figure 1, all NSs are assumed to have the same initial magnetic field strength, 
but different core temperature which do not change as the magnetic field decays. This 
suggests immediately that the key distinguishing property between magnetars and nor- 
mal pulsars is that magnetars were born much hotter than normal pulsars, and thus have 
much longer magnetic field decay time scales, resulting in higher surface magnetic field 
strength even with the same ages of normal pulsars. The above conclusion agrees well 
with the observed correlations between the surface temperatures of magnetars and other 
young NSs, which do not agree with the cooling dominated evolution of neutron stars. 

In the context of temperature evolution and magnetic field decay of magnetized 
NSs, several sophisticated models have developed in recent years (e.g. iGeppert et al. 



20061 ; IPons & Geppertll2007l ; lAguilera et al.ll2008al) . In these studies an isothermal core 
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and a magnetize envelope were generally assumed. They showed that energy released 
by magnetic field decay and Joule heating in the crust is important for the thermal 
evolution of a NS with field strength > lO'^ G. Though the surface temperature of 
radio-quiet pulsars and normal radio pulsars may be broadly covered by the cooling 
curves of different strengthes of initial magnetic field, the surface temperature of mag- 
netars are still too high to reach even the initial field strength is as strong as ~ 10'^ G 
(lAguilera et al.ll2008bL hereafter APM08). 



As pointed in our companion paper dZhang & XidllOlT ). the characteristic ages 



of NSs are typically much longer than the ages of their associated supernova remnants; 
the magnetic field decay is found to be responsible for this apparent discrepancy. Using 
the true ages of NSs, instead of the characteristic ages, will result in reduced core tem- 
peratures for them, by roughly one order of magnitudes; however all other conclusions 
reached above remain unchanged. 

The locations of all magnetars on the left side of the 'death line' in Figure 1, 
in the same way as all other normal pulsars, would suggest that their activities are at 
least related to the mechanism responsible the radiation of normal pulsars. This is 
not expected in the currently accepted model for magnetars, since when their surface 
magnetic fields B exceeds O.lBcr. their pulsar winds m ay be dominated by bou nd pairs 



rather than by freely streaming electrons and positrons (lUsov & Melrosdll996r) . 
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